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1.  INTRODUCTION 


In  Part  I  [Mudaliar,  1993]^  we  obtained  the  dyadic  Green's 
function  (DGF)  of  a  biaxially  anisotropic  medium  for  the  unbounded  case. 
It  may  be  noted  that  the  results  obtained  there  are  in  a  form  very 
convenient  for  deriving  the  DGF's  for  anisotropic  layered  structures. 
Indeed  that  was  one  of  the  underlying  motives  of  Part  I.  We  consider  in 
this  report  the  two -layer  case. 

The  plan  of  the  report  is  as  follows.  In  Section  2  we  briefly 
describe  the  geometry  of  the  problem.  We  formulate  the  solutions  of  the 
problem  in  Section  3.  In  order  to  facilitate  the  evaluation  of  the 
coefficients  of  the  DGF's  we  introduce  amplitude  vectors  and  adopt  a 
matrix  method  in  Section  4.  The  procedure  is  completed  in  Section  5  by 
deriving  various  half- space  Fresnel  coefficients.  Finally  the  report 
concludes  in  Secion  6. 

2.  GEOMMETRY  OF  THE  PROBLEM 


The  geometry  of  the  problem  is  shown  in  Figure  1.  It  consists  of  three 
regions:  Region  0  (  z  >  0  )  is  an  isotropic  medium  with  permittivity 
Region  1  (0>z>— d)  is  the  biaxially  anisotropic  medium  with 
permittivity  c  and  Region  2(z<— d)isan  isotropic  medium  of 
permittivity  c^.  All  the  three  regions  have  the  same  permeability  p.  In 
this  coordinate  system  the  permittivity  c  is  represented  by  the  following 
symmetric  matrix. 
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where  the  elements  c 


ij 


are  explained  In  (3)  of  Part  1. 


3.  FORMULATION 


(1) 


We  have  a  unit  impulse  electric  current  source  in  Region  0  away  from 

the  boundary.  Our  interest  is  to  find  the  DGF's:  and 

GjoCtif').  Here  as  usual  the  second  subscript  denotes  the  region  where 

the  source  is  located  while  the  first  subscript  denotes  the  region 

containing  the  observation  point.  In  order  to  obtain  the  above  DGF's  we 

need  to  solve  the  system  of  equations  given  in  Appendix  A.  Based  on  the 

knowledge  about  the  structure  of  the  DGF  for  the  unbounded  case  obtained 

2 

in  Part  I  and  that  of  the  DGF's  for  the  layered  isotropic  media  we  can 
now  construct  solutions  for  our  two-layer  problem  as  follows. 
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where  for  n  - 


A±  A+ 

and  a  ,b  , 


Che  horizontal  and  vertical  polarizations,  which  are  two  characteristic 

A  A 

linear  polarizations  in  an  isotropic  medium,  a  and  b  denote  the 
polarizations  of  the  a-wave  and  the  b-wave,  which  are  two  characteristic 
(extraordinary)  waves  in  a  biaxially  anisotropic  medium.  Note  Chat  the 
superscripts  +  and  -  indicate  the  upward  and  downward  propagating  waves, 
respectively.  The  task  here  is  to  determine  the  sixteen  reflection  and 
transmission  coefficients:  R's,  A's,  B's  and  X's. 

4.  MATRIX  METHOD 


The  unknown  coefficients  appearing  in  (2),  (3)  and  (4)  can  be 
evaluated  by  using  (A4)  and  (A5) .  But  this  is  at  best  a  tedius  procedure; 
moreover,  the  results  thus  obtained  are  in  a  very  complicated  form.  We 
therefore  follow  an  alternative  method  in  which  we  first  express  the 
unknowns  in  terms  of  half -space  Fresnel  coefficients.  Later  in  the  next 
section  we  evaluate  the  various  half-space  Fresnel  coefficients. 

In  this  procedure  we  denote  the  amplitude  vectors  of  various  waves  as 
p,  q,  P,  Q,  s  (see  Figure  2).  These  are  two-element  vectors  whose 
elements  are  the  characteristic  components  of  the  corresponding  wave.  For 
example 


where  p^  and  p^  are  the  amplitudes  of  horizontally  and  vertically 
polarized  components  of  the  p-wave  in  an  isotropic  medium.  Similarly 
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where  and  are  the  a-  and  b-  conponents  of  the  P-wave  in  a  blaxially 
anisotropic  mediun  and  so  on.  Accordingly  these  amplitude  vectors  satisfy 
the  following  equations 


q 

■^01 

XI 

o 

p 

.^01 

XI 

o 

Q 

■^12' 

- 

P 

s 

.  ^12  . 

(8) 


(9) 


where  ^qj^>  ®-io’  ^10’  ^12  ^12  various  half- space 

reflection  and  transmission  matrices  defined  as  follows 
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*^ab 
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(lib) 


[  ] 


nl2  r  ^  I  *<1.1 

*^aa  **P  L  +k 
=12 

*^ab  ®*P 


’‘av  “P  [  ] 


d12 

Sa 


12 

\h 


[  l(-k“.x*“)d  ]  ■ 
[  ] , 


(12a) 


exp  [  i(-k^^-k2jd  ]  _ 


(12b) 


Note  that  the  exponential  terms  are  included  in  (12)  to  take  into  account 
the  phase  shift  at  the  boundary,  z  -  — d.  For  a  given  incident  wave 
corresponding  to  p,  Equations  (8)  and  (9)  represent  four  equations  for 
four  unknowns  q,  P,  Q  and  s.  The  solution  is  readily  obtained  and  is  given 
as  follows. 


q  -  R  P 


(13a) 
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P  -  A  p 

Q  -  B  p 

s  -  X  p 

where 


(13b) 

(13c) 

(13d) 


(14a) 


(14b) 


(14c) 


(14d) 


We  note  that  through  (14)  we  have  represented  all  the  unknown  coefficients 

in  the  two-layer  DGF's  in  terms  of  half-space  Fresnel  coefficients.  In  the 

uniaxial  limit,  i.e.  when  c  —  c  and  0,  our  results  agree  with  those 

X  y  L 
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of  Lee  and  Kong’.  We  shall  take  up  the  task  of  deriving  explicit 
expressions  for  these  half- space  Fresnel  coefficients  in  the  next  section. 

5.  HALF- SPACE  FRESNEL  COEFFICIENTS 

We  notice  that  for  our  two- layer  problem  there  are  three  situations 
to  consider. 

SITUATION  1 


Here  we  have  free  space  in  the  region  above  the  interface  z  -  0  and 
the  biaxially  anisotropic  medlvun  in  the  region  below  (Figure  2).  A  plane 
wave  from  above  is  incident  on  the  interface.  Depending  on  the 
polarization  of  the  incident  wave  there  are  two  cases  to  consider. 

(1)  h-wave  incidence 


For  this  case  the  electric  fields  in  the  two  regions  may  be 
formulated  as  follows. 


Eo(r)  -  h^  e 


ikQ.f 


01  V  ik  -r 


z  >  0 


(15a) 


Ej^(r) 


.— a  _  ,_b  _ 

A_  Ik  •r  A_  IK  ^r 

a  .  +  b  e 


z  <  0 


The  associated  boundary  conditions  are  given  as  follows . 


(15b) 
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z  X  E^Cr)  -  2  X  Ej^(r)  at  z  -  0  (16a) 

z  X  V  X  1^(2)  -  z  X  V  X  Ej^(r)  at  z  -  0  (16b) 

On  substituting  (15)  in  (16)  we  obtain  the  following  solutions  for  the 


Fresnel  coefficients. 

“to  -  ■  1  \b  ''b- 

C  ^’S.b'V  (17b) 

-  2kor  (  •'x  V  -  ky  V  )  '  (17c) 

*hb  -  ■  7ko3  (  h‘  V  -  h;  yl_  ]  /  i*  (17d) 


where  all  the  unknown  quantities  are  defined  in  Appendix  C. 


(ii)  v-wave  incidence 


Here  the  electric  fields  in  the  two  regions  are  formulated  as  follows. 


_  _  A_  i/c^T  01  A+  ik_T 

E„(r)-Voe  0  .  vj  e  « 


„01 

*  ®vh  "o  " 


z  >  0 


(18a) 
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Ej(r)  - 


va 


A_ 

a 


lie  .r 


A_ 

X  V  b 

vb 


,_b  _ 

ia  ‘r 


z  <  0 


(18b) 


By  subscltuclng  (18)  In  (16),  we  obtain  the  following  solutions. 


w 


- 1  +  X 


va 


(I9a) 


+  X 


vb 


(19b) 


(  Vx  -  Vx  )  C  -  t  Vy  -  Vx  )  V  ]  / 

\b  -  -  [  (  Vx  ■  Vx  )  V  -  (  Vy  -  Vx  )  V  ]  /  ''' 


SITUATION  2 


Here  we  have  the  blaxially  anisotropic  medium  above  the  boundary  z  -  0 
and  the  isotropic  medium  2  below  it.  A  plane  wave  from  above  impinges  on 
the  Interface.  There  are  two  cases  to  consider  depending  upon  the  type  of 
the  incident  wave  . 


(i)  a -wave  incidence 
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The  electric  fields  and  E2  In  the  two  regions  are  formulated 

as 

follows . 

_  A_  lit*«r  12  A+  Ik  ^r  12  Aj.  Ik  ^r 

Ej(r)  -  .  .  +  K..  •  •  *,b  “ 

(20a) 

_  12  ^  1*0  12  ^ 

Ej(r)  -  hj  .  2  -  vj  .  2  .  .  <  0 

(20b) 

Substituting  (20)  In  (16)  we  obtain  the  following  solutions. 

(  V  yb+  -  y>  )  /  i'’ 

(21a) 

“ib  -  •  (  V  xl-  ■  *1-  V  )  ' 

(21b) 

-  R^^v“+  +  R^?v“+  +  v”- 
ah  aa  a^  ab  b^  a 

(21c) 

^  *  V 

- (21d) 

where 

-  k^,  T  1  +  v^,  k.  \A  -  T  , 

fl  yi^  fl  2Z  i  fl^ 

(22) 

for  ?-Cx,y},  f-(a,b},Cl,jl-  C+,-) 

All  other  quantities  are  given  In  Appendix  C;  but  we  have  to  replace 
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(i)  b-wave  incidence 

The  electric  fields  and  E2  In  the  two  regions  are  formulated  as 
follows . 


A_ 

(r)  -  b 


l<c^»r 


12  ^ 


lk^«r 


^  i>12  N- 


lk*»r 


2  >  0 


(23a) 


—  12  ^ 

E,<r)  -  xic  h. 


lK-*r 


.12 


A—  iK-«r 


From  (23)  and  (16)  we  obtain  the  following  solutions. 


Replacements 
a  ^  b 
b  -»  a 


SITUATION  3 


(23b) 


(24) 


In  this  situation  the  geometry  is  the  same  as  in  Situation  1;  the 
only  difference  is  that  the  wave  is  incident  from  below.  Once  again 
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th«re  are  two  cases  to  consider. 


(1)  a-wave  Incidence 


The  electric  fields  In  the  two  regions  are  formulated  as 


“  ^ah  ^0  « 


A+  lk--r  A+  lk-»r 

’  0  4.  Y  «  «  0 


+  X  V.  e 
av  0 


,  z  >  0 


<25a) 


.  -r*  -  b  _ 

—  _  /v^.lk»r  A_lic«r  A_l»c»r 

Ej^(r)  -  a  e  +  a  e  +  be  ,  z  <  0  (25b) 


Applying  the  boundary  conditions  (16),  we  have  the  following  solutions. 


R  -  R^^ 

aa  aa 


12 

ab 

12 

ah 

X  -  X^2 
av  av 


R  ,  -  R 

ab 


X  ,  -  X 

ah 


Replacement  of 
Superscripts 
+  -»  - 

-  -♦  + 


(26) 


Besides  the  above  replacements,  also  note  the  following  changes; 
here  defined  by  (A2)  and  in  other  equations  ^2z  replaced  by  k 


Oz  ■ 


i 


is 


(ii)  b-wave  incidence 
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The  electric  fields  In  the  two  regions  are  formulated  as 


-  *bh  « 


.A+ 


N- 

jv  ^0 


z  >  0 


(27a) 


-  ^  1^‘r  A_  li^.f  A-.  lie^.f 

Ej^(r)  -  be  +R^^be  +R^^ae  .  z  <  0  (27b) 


From  (27)  and  (16)  we  have  the  following  solutions. 

“  *^aa 

^a  "  ^ab 

'S,h  -  ^ah 

*bv  -  V 

It  Is  to  be  noted  that  In  the  uniaxial  limit,  our  results  (17),  (19), 
(21),  (24),  (26)  and  (28)  agree  with  corresponding  results  of  Lee  and 
Kong*. 

6.  NUMERICAL  EXAMPLE 

One  quantity  of  Interest  In  our  results  Is  the  reflection  coefficient 
of  the  two  layer  medium.  To  illustrate  its  characteristics  we  have  taken 
the  following  example.  Region  0  is  free  space  with  permittivity  while 


Replacements 
a  -*  b 
b  a 


(28) 


16 


Region  2  has  pemlttivlcy  (6.0-(-1.006)Cq.  The  anisotropic  medium  is 
characterised  by  the  following  permittivities:  -  (2.8+i.001)cQ  ,  - 

(3.0+i.002)£Q,  -  (3.2+i.003)£Q.  The  tilt  angles  are  -  20*  and  - 

40*.  The  azimuthal  angle  of  incident  wave  is  45*.  The  thickness  of  the 
layer  is  1.0  m  and  the  frequency  is  chosen  as  10  GHz.  In  Figure  3a  we  have 
plotted  the  reflectivities,  r^^  ■  1*^1’  +  ^v  "  ^  ^*^vh^*‘ 

The  behaviour  here  appears  to  be  a  familiar  one.  For  small  angles  on 
incidence  we  notice  that  r^  -  r^  and  there  is  a  null  (Brewster  angle!)  at 
around  60*  for  r^.  The  anisotropic  nature  of  the  medium  is  not  very 
evident  by  looking  at  this  reflectivity  plot.  This  is  perhaps  the  axis  of 
symmetry  is  very  close  to  the  incident  plane.  Next  we  consider  the  case 
when  "  80*  .  Since  the  axis  of  symmetry  is  well  away  from  the  incident 
plane  the  anisotropic  behaviour  is  very  apparent  in  Figure  3b.  Just  to 
illustrate  the  contributions  of  like -polarised  and  cross -polarised 
reflectivities  to  r^  we  have  plotted  )R^|®  and  Figure  3c. 

Perhaps  in  this  example  the  cross -polarised  reflectivity  is  rather  small 
compared  to  the  like  polarised  reflectilvity .  But  the  fact  that  there 
exists  cross-polarization  is  significant. 


7.  CONCLUSION 


We  have  derived  in  this  report  the  dyadic  Green's  function  of  a  two 
layer  biaxially  anisotropic  medium.  The  formulation  is  based  on  the 
solution  for  the  unbounded  case  obtained  in  Part  I .  But  the  evaluation  of 
the  various  coefficients  required  the  use  of  a  matrix  method  which 
expresses  these  two-layer  coefficients  in  terms  of  half-space  reflection 
and  transmission  coefficients.  This  procedure  is  completed  by  deriving 
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explicit  expressions  for  the  various  half-space  Fresnel  reflection  and 
transmission  coefficients.  To  illustrate  the  computational  simplicity  we 
have  provided  a  numerical  example  where  we  have  plotted  the  reflectivities 
versus  the  incident  angle. 
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APPENDIX  A 


The  DGF's  satisfy  the  following  equations 

V  X  V  X  Go„(r.r')  -  J,,  '  i 

7  X  V  X  5jo(r,f')  -  7  Oio(r.f')  -  0  (A2) 

V  X  7  X  Ejoff.r')  -  £j  GjoCr.f)  -  0  (A3) 

where  u  is  the  angular  frequency.  The  boundary  conditions  associated 
with  these  DGF's  are  given  as  follows: 

2  X  GQQ(r,f')  -  2  X  Gj^Q(f,f')  at  2  -  0  (A4a) 

x7xgoo«.f')  -  I  X  7  X  SjgCf.f')  atz-0  (A4b) 

SxGjo(r,f)  -  5xGj„(r,f->  at  z  -  -d  (A5a) 

x7xGjg(r,t-)  -  J  X  7  X  Gj„(r.r')  at  r  -  -d  (A5b) 
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APPENDIX  B 


A+  (  am "  -  A- 

a  ~  [  u  }  X  •  a 


-1/2  -  A+ 


A- 

a 


f  od-^  - 1/2  -  A- 

-  I  J  X  •  <1 


6*  .  ( j .  r 


r  .  ( j .  r 


tu  ^  -2  ^ 

1/  -  t  .  t  , 


t  -  a  or  A 


td  ^  -2  ^ 

1/  «•  t  •  •  t 


t  -  a  or  A 


where 


A+ 

a 


1 

au 


,  a  A 
k  X  o, 


,3.  A 
k  X  o. 


'a  A 


a  A 


|k  X  Oj^l  Ik  X  O2I 


(Bla) 

(Bib) 

(B2a) 

(B2b) 

(B3a) 

(B3b) 

(B4a) 


bu 


C'b  .f'b  A  , 
k  x(k  X  Oj^) 

A  , 
Ik  X  o. 1 


{^b  ,f^b  A  . 
k  x(k  X  o^) 

,^b  A  , 
|k  X  02I 


k^  -  p  k  +  z  k^“ 
P  z 


a  or  b 


(B4b) 

(B5) 


f  (£^xsj.[i^xsj  p/2 

■  1  |kf  X  8, 1  iCf  X  a ,  / 


|k  X  02* 


(B6) 
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A 

k 


f 


-kf 


/  k' 


(B7) 


-  X  (  ±  cos  ^2  +  &2  ^ 


+  y  (  ^  8]^  ^2  ^  ®2  ^1  ^2  ^  ^ 


(B8) 


H  - 


1/2 


(B9a) 


&2 


c 

X 


c 

y 


) 

) 


1/2 


(B9b> 


For  definiteness  we  have  assumed  in  the  above  equations  that  c  <  c  <  t 

X  y  z 


A- 

a 


A 

<c  X  o. 


i^a  A  , 
|k  X  Oj^l 


"a  A 
a  X  o. 


1^*  A  , 

|k  X  <>2* 


(BlOa) 


r 


'^b  /b  A  ^ 
K  X(k  X  Oj^) 

.'^b  A  . 

K  X  O, 


^b  ^'^b  A 
K  X(k  X  O, 


.-^b  A 
a  X  o. 


(BlOb) 


“f  1  ^  A  ,  fd 

a  -  p  k  +  z  k 
P  z 


f  -  a  or  b 


(BID 


yS^  -  {  u  -  d  } 


(B12) 
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APPENDIX  C 


J  4.  f  J  ^  J  i  A  J 
Vi  -  ±  [  »?  r  -  »7  f  J  /  « 

»  X  y  y  X 


where 


{  i.  J  }  -  {  +. 


»»  -  \  h.  V/ 


/•  ^  w 

C  -  \  a,  b  / 


/  h  -♦  + 

»?  -  ') 

I  V  -»  - 


-  .  T  ,  +  .  (-  k-  )  - 

fi  vi$  n  Oz''  £ 


where  ^  |  x,  y} 


,i  i  ,i  i 

0  —  h  V  —  h  V 

X  y  y  X 


T  .  _  k. 

C  'z  z 
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Figure  1 .  Geometry  of  the  Problem 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab; 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology ,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


